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Abstract. The influence of an external field acting differently on the two constituents 
of a binary colloidal mixture performing Brownian dynamics is investigated by 
computer simulations and a simple theory. In our model, one half of the particles 
(A-particles) are pulled by an external force while the other half of them (B- 

particles) are pulled by an external force F^ B \ If and F^ are parallel and the 
field-free state is a mixed fluid, previous simulations (J. Dzubiella et al, Phys. Rev. 
E 65 021402 (2002)) have shown a nonequilibrium pattern formation involving lanes 
of A or B particles only which are sliding against each other in the direction of the 
external forces. In this paper, we generalize the situation both to non-parallel external 
forces and to field-free crystalline states. For non-parallel forces, lane formation is also 
observed but with an orientation tilted with respect to the external forces. If the field- 
free state is crystalline, a continuous increase of the parallel external forces yields a 
novel reentrant freezing behavior: the crystal first melts mechanically via the external 
force and then recrystallizes into demixed crystalline lanes sliding against each other. 
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1. Introduction 

When brought into nonequilibrium, physical systems may spontaneously exhibit many 
different kinds of pattern formation (for recent reviews see Refs. |T], 0) which are much 
richer than the traditional phase transitions in equilibrium systems. While the latter 
are by now well-understood by microscopic theories and simulations || H H, HI, ^ un 
microscopic theories operating on a particle-resolved level for nonequilibrium situations 
still represent a major challenge. In this paper we study a "microscopic" system designed 
to model binary colloidal suspensions in an external field which is an off-lattice version [0] 
of a diffusive system in an external driving field. Colloidal dispersions indeed represent 
excellent model systems which can be brought into nonequilibrium || [| via external 
fields in a controlled way Jl0|| and can be studied experimentally with a real-space 



resolution of an interparticle distance |TT|. Another complementary realization of a 



driven diffusive off-lattice systems is pedestrian motion in pedestrian zones [12, 13, 14] . 

In our model, we consider an equimolar binary mixture of so-called A and B 
particles. The particles are identical as far as their mutual interaction is concerned 
which we model via a Yukawa pair potential having charged colloidal dispersions in 
mind. The A and B particles, however, respond differently to the external field applied: 
A particles feel an external force while I?-particles are driven by a different force 
F^ B \ Completely overdamped Brownian dynamics (with hydrodynamic interactions 
neglected) are assumed |T5| , [T(| for the colloidal motion. The case of parallel forces 
where = —F^ was investigated recently by computer simulation by us and G. 
Hoffmann W7\ via nonequilibrium Brownian dynamics computer simulations in the case 
where the field-free thermodynamic equilibrium state was a mixed fluid. As a result, 
above a critical strength of the external force, the system exhibits a transition towards 
lane formation. The lanes comprise bundles of particles of the same kind (A or B) and 
are parallel to the driving field. This nonequilibrium phase separation ]n| results from 
a slip-stream effect caused or transported by the interparticle interactions. The critical 
force can be theoretically estimated by setting the external force to a typical interaction 
force resulting from a pairwise potential V(r) between the particles. Experimental 
evidence for such an instability has been accumulated in sedimentation dynamics of 
bidisperse suspensions [O, EH, |2T 



In this paper, we generalize the set-up previously studied in Ref. [17] into the 



following two directions: First, we study the case where the two external forces F^ 
and F^ are not parallel. Second, we study - for parallel forces - the case where the 
equilibrium field-free state is crystalline. For non-parallel forces, lane formation does 
also occur. The direction of the lanes, however, is tilted with respect to the driving 
fields. In more details lanes are directed along the difference vector F^ — F^ A ' of 
the two external forces. For a randomly occupied crystal, on the other hand, a two- 
stage transition shows up: first, above a threshold, the external fields melts the solid 
mechanically. Upon increasing the external field strength further, a reentrant freezing 
transition is discovered. The resulting crystalline structure involves completely demixed 
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A- and B solids sliding against each other similar to the fluid lane formation. 

The paper is organized as follows: In section |2], we define the model used and 
describe briefly our simulation technique. Results for non-parallel external forces are 
presented in section |3|. The solid-fluid-solid reentrant behavior is discussed in section |. 
Conclusion are given in section [5[ In particular, we discuss a possible verification of our 
predictions in experiments. 



2. The model 



In our model \T%, we consider a symmetric binary colloidal mixture comprising 2N 



Brownian colloidal particles in d = 2 spatial dimensions. Half of them are particles of 
type A, the other half is of type B with partial number densities Pa = Pb = p/2. The 
colloidal suspension is held at fixed temperature T via the bath of microscopic solvent 
particles. Two colloidal particles are interacting via an effective Yukawa pair potential 

V(r) = Vq aexp [— n(r — o~)/o] /r, (1) 

where r is the center-to-center separation, Vq is an energy scale and a is the particle 
diameter as a length scale. This is a valid model for charge-stabilized suspensions 



confined to two dimensions [£3], ^4|, [2! 

The dynamics of the colloids is overdamped Brownian motion. The friction constant 
£ = 37ir]a (with 77 denoting the shear viscosity of the solvent) is assumed to be the same 
for both A and B particles. The constant external force acting on the zth particle, Fj, is 
different for the both constituents of the binary mixture. It is = F^ for A particles 
and Fi = F^ for B particles. 

The stochastic Langevin equations for the colloidal trajectories fi(t) (i = 1, 2N) 
read as 

^ = E y (\ ^ - ^ i) + p * + ^ (R) (*)■ ( 2 ) 

The right-hand-side includes all forces acting onto the colloidal particles, namely the 
force resulting from inter-particle interactions, the external constant force, and the 
random forces F> describing the collisions of the solvent molecules with the ith colloidal 

particle. The latter are Gaussian random numbers with zero mean, F^ = 0, and 
variance 



(F^) a (t)(Fj K) ) p (t') = 2k B Ti5 aP 8 ij 8{t - t'). (3) 

The subscripts a and f3 stand for the two Cartesian components and ksT is the thermal 
energy. In the absence of an external field, the model reduces to a two-dimensional 
Brownian Yukawa fluid in equilibrium which has been extensively investigated as far as 
structural and dynamical equilibrium correlations and freezing transitions are concerned 



24, 25, 26 



We solve the Langevin equations of motion by Brownian dynamics simulations 
151 ) [2~7| , [28| using a finite time-step and the technique of Ermak 0, p9| . We observed the 
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system running into a steady-state after a transient relaxation. We put N = 250 A and 
iV = 250 B particles into a square cell of length i with periodic boundary conditions. 
The total colloidal number density is p = 2N/£ 2 . The typical size of the time-step was 
0.003r B , where t b = £a 2 /Vo is a suitable Brownian timescale. We simulated typically 
2 x 10 4 time steps which corresponds to a simulation time of 60t#. After an initial 
relaxation period of 20r B , statistics was gathered in the steady state. 



3. Non-parallel external forces 

3.1. General argument for tilted lane formation 

Let us first recapitulate what is know for parallel external forces F^ A ' and : in Ref . 



17J , it was shown that lane formation occurs involving either Aor B particles which are 
sliding against each other in the field direction. In the lane involving A particles only, 
all these particle are drifting with an global velocity while opposite regions which 

involve B particles are streaming with the overall velocity F^ /£. By subtracting the 
overall velocity using a Galilei transformation one readily sees that within the completely 
separated lanes, equilibrium Boltzmann statistics is realized. The system just separates 
into two different equilibrium states which are drifting relative to each other. Physically, 
the formation of lanes is generated by collisions of A against B particles pushed by the 
external force which dynamically separates A and B particles until completely demixed 
lanes are formed. A similar lane formation for sheared granular matter was found via 



molecular dynamics simulation in [30]. The formation of lanes is a sharp first-order 
non-equilibrium phase transition occurring if the external field difference \F^ — F^] 
exceeds a critical value. A Galilei transformation also proves that only the relative 
velocity of A and B regions is relevant. Hence, without loss of generality, it is sufficient 
to study the special case F^ A ' = —F^ B '. 

For non-parallel external forces, the collisions between A and B particles are not any 
longer central and the phase separated structure will be different, in general. In order 
to get insight into the location of an interface between two completely demixed regions 
involving A and B particles only, we first do a simple continuum argument: consider a 
(one-dimensional) interface between an A and B region with a direction described by a 
two-dimensional vector d, see Fig. |l|. The full interface position can be parameterized 
by a set of vectors 



f(s) = R + sd, (4) 

where s is a real parameter and Rq is a vector describing a point on the interface. After 
a time At, the A particles have moved on average a distance F^At/^, while the B 
particles have been displaced by a distance F^At/^. Neglecting any collisions, the A 
particles near the interface will move towards a new interface which is described by the 
set of vectors 

r< A) (s) =R + sd + F {A) At/^ (5) 
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Figure 1: One-dimensional interface separating a region containing A particles (light 
gray) and B particles only (dark gray). The direction of the interface is d. The vector 
Ro points to the interface. The new interface after a time At at distance V\At from the 
original one is shown ClS db dashed line. 



while the B particles near the interface will drift to 

r< B) (s') =R + s'd+ F {B) At/f (6) 

with another real parameter s' . The interface can only be stable if these two interfaces 
coincide. In case they do not, there is either empty space which will be filled by 
neighboring particles or A and B particles will collide which will destroy the interface 
as well. Hence, the stability criterion is that for any s there is an s' such that 
r^ A '{s) = f^ B \s') which simply yields the condition 



d 



At 



This implies that for an interface to be stable, its direction has to be collinear with the 
force difference 

p(B) _ p(A)_ 



AF 

Hence, the angle a describing the interface orientation (see Fig. |1|) is 

AF • F^ 



a 



arcsm • 



(8) 



(9) 



|AF||F( fi )| 

Clearly, contrarily to the case of parallel forces, the interface position will move in space 
The interface velocity vi normal to its position can be calculated as 

\p(A) x p(B)\ 



Vi 



f|AF| 



(10) 
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Obviously, the same argument can be repeated with exchanged roles of A and B showing 
that stable parallel lanes with the direction AF are expected with move with the 
interface velocity vi given by Eq. ([10]). Furthermore, the same argument applied for 
a small At shows that a curved interface is unstable such that a stable interface has to 
be straight. 

Let us finally discuss two special cases: first returning to parallel forces, indeed the 
interface direction is parallel to the field direction, the angle a is tt/2, and the interface 
velocity vanishes as follows directly from Eq. (|10"D . Second the case of perpendicular 
forces deserves some particular attention. Here the angle a is 

a = arctan — (11) 

\FW\ 

and the interface velocity can be expressed as 

vi = - l± ^ J -. 12 

ciafi 

Obviously, this general argument is only a necessary condition for a stable AB 
interface. An alternative is a mixed situation with no interface at all driven by entropy. 
In analogy to the parallel case we anticipate that a critical strength of AF is needed 
to build stable lanes. This will now be investigated by computer simulation and simple 
theory in more detail. 

3.2. Computer simulation results 

In our simulation we assume - without loss of generality - the direction of along 
the y axis of the simulation box. The model parameters are fixed to pa 2 = 1.0, k = 4.0, 
and Vq = 2.5/cbT. Simulation snapshots for different external field strengths are shown 
in Fig. l(a)-(d). 

In Fig. 0(a), no field is applied and a homogeneous completely mixed state is visible. 
In Fig. |H(b), on the other hand, the external forces are parallel: F^ = —F^ B \ The 
magnitude \F^ — F( B '\ is beyond the critical strength JlTj] such that lane formation 



parallel to the field shows up. Nonparallel forces with perpendicular directions are 
investigated in Figs. |2] (c) and (d). One observes formation of tilted lanes which 
are indeed in the direction of the force difference vector as expected from our general 
argument. In following the configurations as a function of time we verified the simple 
formula of the interface velocity vi as given in Eq. flTOl). 

In a next step, we study perpendicular forces keeping their ratio q = \F^\/\F^ B >\ 
fixed but increasing their magnitude. For small forces the system stays demixed. In 
increasing the strength of both fields, we have calculated suitable order parameters 
which are sensitive to tilted lane formation. These are immediate generalizations of 
those used in Ref. [I7J]. The order parameter exhibits a sharp jump indicating the 



critical field strength of the force difference vector upon which tilted lane formation is 
achieved. There is a clear hysteresis loop if the force is reduced again such that the non- 
equilibrium phase transition towards tilted lane formation is of first order. The critical 
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Figure 2: Typical simulation snapshots of the two dimensional system: (a) disordered 
state without field, (b) lane formation with parallel fields F^ A > = — = 180£;bT/<7 
in y-direction above the critical force, (c) lane formation with perpendicular fields of 
same magnitude = = 180/csT/cr above the critical force difference, (d) lane 

formation with perpendicular fields and |F (A) |/I^ (B) I = V 2 with = 180fc B T/a 

above the critical force difference. In (c) and (d) the A particles are drifting in 
x-direction, while B particles are drifting in ^-direction. The lanes are moving 
perpendicular to AF. The particles are depicted as spheres with diameter a. A light 
sphere is an A-particle while a gray sphere is a 5-particle. 

value AF C of the force difference | AF\ is shown versus the ratio q = \F^\/\F^\ in Fig. 
|3|. We shall compare these data to a simple theoretical prediction in the next subsection. 

A careful remark is in order for small positive q. The periodic boundary conditions 
used in the simulation corresponds to a toroidal topology shown in Fig. f|. If q is getting 
small the boundary conditions enforce a multiple winding around the torus such that 
finite size effects are expected to be significant. Therefore we have not shown simulation 
data for small q in Fig. 
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Figure 3: Dimensionless critical force difference AF* = AF c a/kBT versus the ratio 
q between the perpendicular forces and \F^ B ^ for the two-dimensional system. 

Above the critical force the system is in the patterned state characterized by stripe 
formation. The circles are simulation results, while the solid line is theory. The 
parameters are k = 4.0, V = 2.5k B T, and pa 2 = 1.0. 



3.3. Theory 

We are aiming at a rough theoretical estimation of the boundaries of the demixing 
transition with constant external fields F^ and F^ B \ This was already put forward 



in Ref. [17] for parallel forces and is generalized here to the general case of non-parallel 
forces. Consider first a central collision between an A and B particle pair, see Fig. 
|5|. Transforming the trajectories onto one with a fixed common center of mass of the 
two particles, one realizes that the collision is effectively driven by half of the force 
difference, \{F^ - F^) = -|AF for A particles and \{F^ - F^) = |AF for B 
particles. A transition towards patterned lanes is expected if |AF|/2 is larger than a 
typical average force between A and B particles, and lane formation is induced. The 
latter force depends both on density and on the external fields themselves. We estimate 
a typical average force between two opposite particles by considering different "effective" 
interparticle spacings. The first typical interparticle spacing is set by the density alone, 
a = p~ x l 2 . Including fluctuations in the interparticle distance induced by a finite 
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Figure 4: Visualization of the two-dimensional system by mapping a two-dimensional 
simulation snapshot onto the surface of a torus. A particles are drifting along the torus, 
while B particles are drifting around it. The lanes are moving perpendicular to their 
direction. 




Figure 5: Sketch of two colloids of opposite type colliding due to the external fields 
and F( B \ The A particle is white, while the B particle is gray. The particles collide 
effectively in direction of the difference force AF = F^ — F^ as seen from a fixed 
center-of-mass. 
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temperature results in a further smaller effective average distance a as obtained by 
setting a typical interparticle energy equal to V(a) + k B T. Hence a = V~ x [V(a) + ksT] 
where V~ x is the inverse function of the interaction potential V(r). Finally the presence 
of the external fields enforces an even smaller averaged distance a' between two colliding 
opposite particles. 

We estimate this minimum distance a' by adding the net force per colliding particle 
AF/2 to the force at distance a via 

a' = F- 1 [AF/2 + F(a)] , (13) 

where F _1 is the inverse function of F(r) = — ^V( r )- In general, an AB particle pair 
will not collide directly along AF/2 such that the actual average distance is between a' 
and a. Hence the averaged force / between an A and a B particle is roughly 

7 = ^[V(a')-V(a)}. (14) 
a — a 

The critical force difference AF C is reached when it becomes of the order of the mean 
force /, 

AF C = 2XJ. (15) 
A = 2 is a dimensionless prefactor which is determined by an optimal fit to all simulation 



results for parallel forces |L7 



The basic prediction of this simple theory is that the only essential parameter 
governing tilted lane formation is the magnitude |AF|. This prediction can be tested 
by simulation. In fact, in Fig. [| the critical value of \AF\ is shown for different 
ratios q = \F^\/\F^\. Were the theory correct, all the simulation data should fall 
on a horizontal line independent of q. As can be deduced from Fig. |3], this is indeed 
confirmed. Furthermore, the actual magnitude predicted from the theory is in line with 
the simulation data. Note that - as far as the non-parallel case is concerned - there is 
no fit parameter involved. The global fit parameter A is solely adjusted to the case of 
parallel forces. 



S.Jf.. Three-dimensional model 

The model and all methods and arguments can readily be generalized to three spatial 
dimensions. Similar conclusions hold for the formation of tilted lanes. We have also 
performed computer simulations in a cubic box in three dimensions and observed tilted 
lane formation. Results are presented in Fig. [| tilted lane formation is clearly visible in 
the plane spanned by the two forces F^ and F^ B \ see Fig. |6| (a). Perpendicular to the 
direction AF of the lanes, the system shows a structure reminiscent of two-dimensional 
spinodal decomposition, see Fig. || (b). The parameters are for these snapshots k = 4.0, 
Vo = 2.5/c#T and pa 3 = 1. In conclusion, this shows that pattern formation is a general 
effect which is independent of the dimensionality of the model. 




Figure 6: Typical snapshots of the three-dimensional system with perpendicular external 
fields above the critical force difference. The mag nitude of the forces is \F^\ = \F^\ = 
150k bT /a. (a) Three-dimensional view, (b) look on the plane perpendicular to the AF 
vector, cd 



4. Reentrant freezing for a driven Brownian crystal 

In this section we focus on a Brownian crystal which is driven by an external field. To this 
end, the external fields acting onto A and B particles are parallel (F^ = —F^ = fe y ) 



but in contrast to Ref. |17j] the field-free equilibrium initial state is a triangular crystal. 
The case of a fluid field-free configuration is easier as rotational symmetry with respect 
to the direction of the external field applied is ensured. This is no longer true for 
a crystal where one has to specify the field direction with respect to the crystalline 
orientation resulting in an anisotropy which reflects the crystalline symmetry. Clearly, 
in the absence of any external field, A and B particles are indistinguishable; hence the 
equilibrium state is randomly occupied triangular crystal. In the other limit of very 
strong fields, one expects again phase separation into completely demixed A and B 
regions. Once they are demixed, they follow Boltzmann statistics. Consequently the 
equilibrium state is a pure A (or B) crystal of the same lattice than the original (field 
free) one. What is less clear intuitively is how the system transforms from the first 
randomly occupied crystal into the demixed crystal if the field is turned continuously 
on. At least two scenarios are conceivable: either the system retains the underlying solid 
lattice but particle exchange hopping processes generated by the external field demix 
the crystalline state or the crystal first melts mechanically via the external field and 
then crystallizes again. In our simulations we almost exclusively observed the latter 
scenario. 
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Figure 7: Crystallinity order parameter \]/ 6 (solid line) and lane order parameter <ft 
(dashed line) plotted versus reduced external force /* as calculated by a Brownian 
dynamics simulation. 



In order to detect a triangular crystalline order we define a suitable crystallinity 
order parameter that probes sixfold symmetry around a given particle via [03], |3T 

2N 



1 



EE 



,6i</>jfc 



(16) 



\ 12iV^ 

\ 3=1 <k> 

Here the fc-sum includes the six nearest neighbors of the given particle and the j-sum 
extends over 2N particles in the simulation box. The large angular brackets indicate 
a time average, ipjk is the polar angle of the interparticle distance vector with respect 
to a fixed reference frame. For ideal sixfold symmetry, i.e., for a perfect triangular 
crystal, — 1- Thermal fluctuations cause deviations from this ideal case but a value 
of ^6 > 0.8 [[H], [3^] is conveniently taken to be conclusive evidence for a triangular 
crystal. 

In Fig. [7] we plot \l/ 6 versus the difference external force /* = \F^ B >\a /ksT for 
fixed parameters k = 4.0, Vq = lb.OksT, and pa 2 = 1. The direction of the force 
is (ll)-direction of the triangular crystal. Note that in contrast to the parameters 
used in section ^, the interactions energy is much larger to ensure that the equilibrium 
field-free state is crystalline. Indeed the crystallinity order parameter in the field-free 
case is around = 0.87. Upon increasing the external field strength to /* ~ 10, 
the crystallinity order parameter sharply drops down and stays to values close to zero 
indicating a melting of the crystal induced by the external field. This melting process 
is mainly caused by a mechanical stress induced by the external field with respect to 
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Figure 8: Snapshots of the two-dimensional system for different external forces, starting 
with a solid in the free-field state. The forces acting on the two different particle types 
are pointing into opposite directions. The forces are (a) /* = 0, (b) /* = 50, (c) 
/* = 150, (d) /* = 250. In (a) and (d), the system shows a solid structure, while in 
(b) and (c) the system is a homogeneously fluid. Here the total number of particles is 
2N = 250. 



AB particle pairs. As the field is getting larger (/* > 200), the order parameter \&6 
increases again to values close to 0.7. This is accompanied with lane formation as 
indicated by a drastic increase of another order parameter defined in Ref. [17] which 
is sensitive to lane formation. The whole scenario is illustrated also by simulation 
snapshots shown in Fig. |8|. While Fig. || (a) corresponds to a field-free randomly 
mixed crystal, Fig. || (b) and (c) are in the molten state while Fig. [| (d) represents 
a refrozen demixed crystal sliding against each other. In Fig. [8] (b) and (c) worm like 
structures along the fields occur as a precursor to lanes formed by solids. Consequently 
we have shown evidence for a reentrant freezing behavior generated by external fields 
in non-equilibrium. A qualitative similar situation occurs for colloidal solids in linear 
shear flow |33| , |34] , |35] . A continuous increase of the shear rate can lead to shear- melting 
and subsequent recrystallization into a different solid structure |36], [37], [38]]. A similar 



effect is shear thinning and subsequent shear thickening as observed in colloidal fluids 
for increasing shear rates |39| . 
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5. Conclusions 

To summarize: we have generalized earlier studies of non-equilibrium pattern formation 
in continuum driven diffusive mixtures to non-parallel external fields and crystalline 
states. As main results we found tilted lanes along the force difference vector which are 
wandering with a constant interface velocity provided the external force difference is 
large enough. Furthermore, a solid melts and refreezes if the magnitude of an external 
field is increased. 

Let us remark on possibilities of non-equilibrium lane formation in more 
complicated systems: first, our system studied was completely symmetric involving 
the same partial densities and the same particle-particle interactions. For experimental 



realizations ||T9| , p20| , [21|j this will not be fulfilled in general. However, the basic physics 
of lane formation will not change. Secondly, if ternary and further multicomponent 
mixture beyond binary ones are considered, we expect cascades of lane formation 
transitions involving the different particles species as the external field is increased. 

We finally comment on possible experimental realizations of our model: There 
are different fields where the pattern formation we predicted within our model can 
be verified, namely in colloidal dynamics and in pedestrian motion. Binary colloidal 
mixtures indeed can be driven by constant external forces. Important examples 
for parallel forces are sedimentation where the external force is gravity |19|, |20| or 



electrokinetic motion of charged colloids [HDJ where the external force is an electric 



field. Both the fluid and crystalline field-free state can be studied. A recent realization 
with mixed crystals can be found in |41||. One possible drawback is the hydrodynamic 



backflow caused by strong hydrodynamic interactions |43|, (44], [45[ which were 
neglected in our model. An overall backflow can be avoided by a time dependent 
oscillatory field (e.g. AC electric field) which leads qualitatively to the same lane 
formation if its frequency is small enough [lTj. We further think that the long-ranged 
hydrodynamic flow around a driven colloidal particle will favor lane formation, i.e. the 
critical field strength needed to generate lane formation is expected to be lower than 
with hydrodynamic interactions neglected. Colloids can also be exposed to external 
laser-optical and magnetic filed |H| which generate external forces in a controlled way 
coupling to the dielectricity (resp. the magnetic permeability) of the colloidal material. 
Nonparallel external forces in colloidal mixtures can be realized by crossing two external 
fields e.g. gravity with electric, laser-optical with electric, laser-optical with magnetic 
etc. The two species of a colloidal mixture will in general respond differently to the two 
external fields such that the resulting total external force will be different in direction. 

A different realization is dynamics of pedestrians in pedestrian zones and in lecture 
halls. Similar off-lattice models involving Brownian particles have been used to simulate 
the collective behavior of pedestrians [T3| including escape panic f£E| . Our setup 
of perpendicular external fields is realized by two crossing pedestrian lanes in which 
pedestrians are only moving in one direction. Based on our results, we would expect 
tilted lane formation provided the density of the pedestrians is high enough. Finally it 
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would be interesting to extend phenomenological hydro dynamical theories which predict 
lane formation for parallel forces via an instability |57], ^8], [I 
of tilted forces. Work along these lines is in progress. 
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